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We show that it is impossible to determine the time a tunneling particle spends under the 
barrier. However, it is possible to determine the asymptotic time, i.e., the time the particle spends 
in a large area including the barrier. We propose a model of time measurements. The model 
provides a procedure for calculation of the asymptotic tunneling and reflection times. The model 
also demonstrates the impossibility of determination of the time the tunneling particle spends under 
the barrier. Examples for J-form and rectangular barrier illustrate the obtained results. 
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Tunneling phenomena are inherent in numerous quantum systems, ranging from atom and condensed matter to 
quantum fields. Therefore, the questions about the tunneling mechanisms are important. There have been many 
attempts to define a physical time for tunneling processes, since this question has been raised by MacColl [Q in 
1932. This question is still the subject of much controversy, since numerous theories contradict each other in their 
\ predictions for "the tunneling time". Some of these theories predict that the tunneling process is faster than light, 
^^O ■ whereas the others state that it should be subluminal. This subject has been covered in a number of reviews (Hauge 
and St0vneng ||, 1989 ; Olkholovsky and Recami ||, 1992; Landauer and Martin §, 1994 and Chiao and Steinberg 
||, 1997). The fact that there is a time related to the tunneling process has been observed experimentally [p[-^4|. 
However, the results of the experiments are ambiguous. 
, Many of the theoretical approaches can be divided into three categories. First, one can study evolution of the wave 
packets through the barrier and get the phase time. However, the correctness of the definition of this time is highly 
questionable Jig ]. Another approach is based on the determination of a set of dynamic paths, i.e., the calculation of 
the time the different paths spend in the barrier and averaging over the set of the paths. The paths can be found from 
the Feynman path integral formalism jl(|, from the Bohm approach [[l7|-fj0||, or from the Wigner distribution [glj. 
The third class uses a physical clock which can be used for determination of the time elapsed during the tunneling 
(Biittiker and Landauer used an oscillatory barrier [[l5), Baz' suggested the Larmor time p2|). 

The problems rise also from the fact, that the arrival time of a particle to the definite spatial point is a classical 
concept. Its quantum counterpart is problematic, even for the free particle case. In classical mechanics for the 
determination of the time the particle spends moving along a certain trajectory we have to measure the position of 
the particle at two different moments of time. In quantum mechanics this procedure does not work. From Heisenberg's 
uncertainty principle it follows that we cannot measure the position of a particle without alteration of its momentum. 
To determine exactly the arrival time of a particle, one has to measure the position of the particle with a great 
precision. Because of the measurement, the momentum of the particle will have a big uncertainty and the second 
measurement will be indefinite. If we want to ask about the time in quantum mechanics, we need to define the 
procedure of the measurement. We can measure the position of the particle only with a finite precision and get a 
distribution of the possible positions. Applying such a measurement, we can expect to obtain not a single value of 
the traversal time but a distribution of times. 

The question How much time does the tunneling particle spend in the barrier region? is not precise. There are two 
different, but related questions connected with the tunneling time problem 

(i) How much time does the tunneling particle spend under the barrier? 

(ii) At what time does the particle arrive to the point behind the barrier? 

There have been many attempts to answer these questions. However, there are several papers showing, that 
according to quantum mechanics the question (i) makes no sense ]23|-p6| . The goal of this paper is to investigate the 
possibility to determine the tunneling time using a concrete model of time measurements. 

The paper is organized as follows: In secti on [ll| we prove that it is impossible to determine the time the tunneling 
particle spends under the barrier. In section IHTwe present the procedure of the time measurement. This procedure 
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leads to the dwell time if no distinctions between tunneled and reflected particles are made. This is shown in section 
IV . In section [v] we modify the proposed procedure of the time measurement to make the distinction between tunneled 
and reflected particles and obtain the tunneling time. The result of such a procedure clearly shows the impossibility 
of the determin ation of the tunneling time. However, it also gives the method of the asymptotic time calculation. In 
sections VI and VII we examine the pro perties of the tunneling and reflection times. In section 
formula for asymptotic time. Section 



[X summarizes our findings. 



VIII, we derive the 



II. IMPOSSIBILITY OF THE TUNNELING TIME DETERMINATION 



To answer the question How much time does the tunneling particle spend under the barrier? we need a criterion of 
the tunneling. In this paper we accept the following criterion: the particle had tunneled in the case it was in front 
of the barrier at first and later it was found behind the barrier. We require that the mean energy of the particle and 
the energy uncertainty must be less than the height of the barrier. Following this criterion, we introduce an operator 
corresponding to the "tunneling-flag" observable. This operator projects the wave function onto the subspace of 
functions localized behind the barrier 



h (x) = e(x-x) 



(i) 



where is the Heaviside unit step function and X is a point behind the barrier. We call operator fa as the tunneling 
flag operator. This operator has two eigenvalues: and 1. The eigenvalue corresponds to the fact that the particle 
has not tunneled, while the eigenvalue 1 corresponds to the tunneled particle. 

We will work in the Heisenberg's representation. In this representation, the tunneling flag operator is 



fT(t,X) = exp[ ^Ht\ fa(X)exp(-^Ht 



(2) 



To take into account all tunneled particles, the limit t — -> +oo must be taken. So, the "tunneling-flag" observable 
in the Heisenberg's picture is represented by the operator fa (oo,X) = limi^ +00 fa (t,X). We can obtain explicit 
expression for this operator. 

The operator fa (t, X) obeys the equation 



ihj- t fa(t,x)= /, u.xi.u 



(3) 



The commutator in Eq. 



exp I —Ht 

n 



fa(X),H 



exp 



-Ht 



may be expressed as 
fa(t,X),H = 

If the Hamiltonian has the form H = ^p 2 + V (x) then the commutator takes the form 

fa (X) , ff] = ihJ (X) 

where J (X) is the probability flux operator 

^ X ) = 2M ^ ( X \P + P \ x ) ^l) ■ 
Therefore, we have equation for the commutator 

/ T (i,x),£] =ihJ(x,t). 

The initial condition for the function / (t, X) may be defined as 

fa(t = 0,X) = fa(X). 
From Eqs. (0) and (p|) we obtain equation for the evolution of the tunneling flag operator 



(4) 



(5) 



(6) 
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%nj- t f T {t,x) = %nj{x,t). (7) 

From Eq. (Q) and initial condition it follows an explicit expression for the tunneling flag operator 

f T (t,X) = f T (X)+[dt 1 J(X,t 1 ). (8) 
Jo 

In the question "How much time docs the tunneling particle spend under the barrier?" we ask about the particles, 
which wc know with certainty have tunneled. In addition, we want to have some information about the location of 
the particle. However, does quantum mechanics allow us to have the information about the tunneling and location 
simultaneously? A projection operator 

D(T) = J dx\x) (x\ (9) 
r 

where \x) is the eigenfunction of the coordinate operator, represents the probability for the particle to be in the region 
r. In the Heisenberg's representation this operator takes the form 

D (r, t) = exp f^Htj D (r) exp (- ~Ht\ . (10) 

From Eqs. (§), (§ and ((TJ) we see that the operators D(T,t) and /t(oo,A) in general do not commute. This 
means that we cannot simultaneously have the information about the tunneling and location of the particle. If wc 
know with certainty that the particle has tunneled then we can say nothing about its location in the past and if we 
know something about the location of the particle, we cannot determine definitely whether the particle will tunnel. 
Therefore, the question How much time does the tunneling particle spend under the barrier? cannot be answered 
in principle, if the question is so posed that its precise definition requires the existence of the joint probability that 
the particle is found in T at time t and whether or not it is found on the right side of the barrier at a sufficiently 
later time. A similar analysis has been performed in Ref. ( p6| . It has been shown that due to non-commutability of 
operator there exist no unique decomposition of the dwell time. 

This conclusion is, however, not only negative. We know that J_ °° dx \x) (x\ = 1 and 1, fp (oo, X) = 0. Therefore, 

if the region T is large enough, one has a possibility to answer the question about the tunneling time. 

From the fact that the operators D (T, t) and fp (oo, X) do not commute we can predict that the measurement of the 
tunneling time will yield a value dependent on the particular detection scheme. The detector is made so that it yields 
some value. But if we try to measure non-commuting observables, the measured values depends on the interaction 
between the detector and the measured system. So, in the definition of the Larmor time there is a dependence on the 
type of boundary attributed to the magnetic-field region (3) . 



III. THE MODEL OF THE TIME MEASUREMENT 



We consider a model for the tunneling time measurement which is somewhat similar to the "gedanken" experiment 
used to obtain the Larmor time, but it is simpler and more transparent. This model had been proposed by A. M. 
Steinberg however, it was treated in non-standard way, introducing complex probabilities. Here we use only the 
formalism of the standard quantum mechanics. 

Our system consists of particle P and several detectors D. Each detector interacts with the particle only in the 
narrow region of space. The configuration of the system is shown in Fig. [|. When the interaction of the particle with 
the detectors is weak the detectors do not influence the state of the particle. Therefore, we can analyze the action of 
detectors separately. 

First of all we consider the interaction of the particle with one detector. The Hamiltonian of the system is 

H = Hp + Hp, + Hj (11) 
where Hp = -^jp 2 + V (x) is the Hamiltonian of the particle, Hp, is the detector's Hamiltonian and 

H I = 1 qD(x D ) (12) 
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FIG. 1. The configuration of the measurements of the tunneling time. The particle P is tunneling along x coordinate and it 
is interacting with detectors D. The barrier is represented by the rectangle. The interaction with the definite detector occurs 
only in the narrow region limited by the horizontal lines. The changes of the momenta of the detectors are represented by 
arrows. 



represents the interaction between the particle and the detector. The operator q acts in the Hilbcrt space of the 
detector. We require a continuous spectrum of the operator q. For simplicity, we can consider this operator as 
the coordinate of the detector. The operator D (xd) acts in the Hilbert space of the particle. In the coordinate 
representation it is nonvanishing only in the small region around the point xd- In an ideal case the operator D (xd) 
may be expressed as <5-function of the particle coordinate, 

D (x D ) = \x D ) (x D \ = 5 (x - x D ) ■ (13) 

Parameter 7 in Eq. ( p"2| ) characterizes the strength of the interaction. A very small parameter 7 ensures the undis- 
turbance of the particle's motion. 

Hamiltonian ([l2]) with D given by (|l^) represents the constant force acting on the detector D when the particle 
is very close to the point xd- This force results in the change of the detector's momentum. From the classical point 
of view, the change of the momentum is proportional to the time the particle spends in the region around Xd and 
the coefficient of proportionality equals to the force acting on the detector. In the ordinary quantum mechanics there 
is no general method of the time determination. If we want to define such a method, we have to make additional 
assumptions about the time. It is natural to extend the classical method of the time determination into the quantum 
mechanics, too. Therefore we assume that the change of the mean momentum of the detector is proportional to the 
time the constant force acts on the detector and that the time the particle spends in the detector's region is the same 
as the time the force acts on the detector. 

We can replace the (5-function by the narrow rectangle of the width L and height . From (|l2|) it follows that the 
force acting on the detector when the particle is in the region around xd is ~7jj- The time the particle spends in the 

region around xo equals to (— 7r) 1 ((p q (t)) — (Pq)) where p q is the momentum of the detector conjugated to the 
coordinate q, while (p q ) and (p q (t)) are the mean initial momentum and momentum after time t, respectively. The 
time the particle spends until time moment t in the unit length region is 

r(t) = --({p q (t))-{p q )). (14) 
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To find the time the particle spends in the region of the finite length, we have to add the times spent in the regions 
of length L. When L^Owe obtain an integral. 
The evolution operator is 



U (t) = exp 



1 



(Hp + H d + JT/) t 



(15) 



In the moment t = the density matrix of the whole system is p(0) = pp (0) £g> pp> (0) where pp (0) is the density 
matrix of the particle and po (0) = |3>) ($| is the density matrix of the detector with |$) being the normalized 
vector in the Hilbert space of the detector. After the interaction, the density matrix of the detector is pp, (t) — 

Tr P (ft (t) (pp (0) <g> |$) ($|) U + (t)\. In the moment t = it must be (a;| p P (0) |a;') ^ only when x < and x' < 0. 

Further, for simplicity we will neglect the Hamiltonian of the detector. The evolution operator then approximately 
equals to the operator U (t, jq) where 



U (t, a) = exp 



-(— 



p 2 + V(x) + aD (x D ) t 



(16) 



After such assumptions from our model we can obtain the time, the particle spends in the definite space region. 
Similar calculations were done for detector's position rather than momentum by G. Iannaccone p8fl. 



IV. MEASUREMENT OF THE DWELL TIME 



We expand the operator U (t,jq) into the series of the parameter 7, assuming that 7 is small. Introducing the 
operator D (xd) in the interaction representation 

(17) 



(18) 



D (x D ,t) = exp y-Hptj D (x D ) exp ( -^H P t 
we obtain the first order approximation for the operator U (t,jq), 



U{t,jq) nexp(-^H P t) U + ^ J dt x D {x D M) 
For shortening the notation we introduce an operator 

F(x D ,t)= f dtxb{x D ,t x ) 
Jo 

and the equation for the evolution operator U (i, jq) is expressed as 

79 



t/(t,7g)«exp^--# P tJ ^l + J^F(x D ,t) 
The density matrix of the detector in the coordinate representation in the first order approximation then is 
(q\ PD (t) W) = (q m <$ W) Tr [ll (t, jq) p P (0) & (t, W)} 

= (q |$) <$ \q') (l + n(p{ XD ,t))-^ (P (x D ,t) 



(19) 



(20) 



(q\ exp 



j{F(x D ,t))q 



|$) ($| exp 



1 (F{x Dl t))q 



The average momentum of the detector after time t is (p q ) — 7 (^F (xo,t)j where ( p q ) = ($1 p q 1$) and (^F (xd, t) 

Tr ^F (xd, t) pp (0)|. From Eq. ( p^ ) we obtain the time the particle spends in the unit length region between time 
momentum t = and t 

(21) 



r Dw {x,t) = (F{x,t) 
The time spent in the space region restricted by the coordinates x\ and X2 is 



t Dw (x 2 ,x 1 )= dxr Bw (x,t-> 00) = 



dx 



p (x, t) dt. 



(22) 



This is a well-known expression for the dwell time || . The dwell time is the average over entire ensemble of particles 
regardless they are tunneled or not. The expression for the dwell time obtained in our model is the same as the 
well-known expression obtained by other authors. Therefore, we can expect that our model can yield a reasonable 
expression for the tunneling time as well. 
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V. CONDITIONAL PROBABILITIES AND THE TUNNELING TIME 



Having seen that our model gives the time averaged over the entire ensemble of particles, let us now take the average 
only over the subensemble of the tunneled particles. The joint probability that the particle has tunneled and the 

detector has the momentum p q at the time moment t is W (fT,Pq]t) — Tr |/t (X) \p q ) (p q \ p(i)| where \p q ) is the 

eigenfunction of the momentum operator p q and the tunneling flag operator fx [X] is defined by Eq. (|l]). In quantum 
mechanics such a probability does not always exist. If the joint probability does not exist then the concept of the 
conditional probability is meaningless. But in our case the operators Jt (X) and \p q ) (p q \ commute, therefore, the 
probability W (fT,p q ',t) exists. The conditional probability that the momentum of the detector is p q provided that 
the particle has tunneled is given according to the Bayes's theorem, i.e., 



W{p q ;t\fr 



W(f T ,p q ;t) 
W(f T ;t) 



(23) 



where W (fr',t) — Tr |/t (X) p (i)| is the probability that the particle has tunneled until time t. The average 
momentum of the detector with the condition that the particle has tunneled is 

(p,(t)) = / Pq d Pq W (p,;*|/ T ) 



or 



1 



-Tr 



{f T {X)p q p{t)). 



W(f T ; t) 

In the first order approximation the probability W (/x; t) is given by equation 

W (/t; t) » (j T (t, X)) + ^ (q) ( [f T (t, X) , F (x D ,t) 

The expression Tr |/t (X)p q p (i) j in Eq. (p^)in the first order approximation reads 

Tr[f T (X) Pq p(t)} «<p g )(/ T (*,X)J 



(24) 



(25) 



+ ^ ((It (t, X) F (xn,t)) (p q q) - (q Pq ) (f (x D ,t) f T (t, X) 



Using the commutator [<7,p 9 ] = %h from Eqs. ( |l4] ) and (|2j) we obtain the time the tunneled particle spends in the 
unit length region around x until time t 



r(x,t) 



1 



2{f T (t,X] 
1 



ih(f T (t,x) 



f T (t,X)F(x,t) + F(x,t)f T {t,X) 
((q) (p q ) - Re (qp q ) ) ( \f T (t, X) , F {X, t) 



(26) 



The obtained expression ( ^6] ) for the tunneling time is real, on the contrary to the complex-time approach. It should 
be noted that this expression even in the limit of the very weak measurement depends on the particular detector. 

This yields from the impossibility of the determination of the tunneling time. If the commutator fx {t, X) , F (x, t) 

is zero, the time has a precise value. If the commutator is not zero, only the integral of this e xpres sion over a large 
region has the meaning of an asymptotic time related to the large region as we will see in Sec. VIII . 

Eq. ( |26] ) can be rewritten as a sum of two terms, the first term being independent of the detector and the second 
dependent, i.e., 



r(x,t) 



„Tun 



(x, t) + - ((<?) {p q ) - Re (q Pq )) t c ^ (x, t) 



(27) 



where 
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T corr i X l ^) — 



1 



1 

2i(f T (t,X) 



f T (t,X)F(x,t) + F(x,t)f T (t,X) 
f T {t,X),F(x,t) 



(28a) 
(28b) 



The quantities r^ 11 (x,i) and (x, t) do not depend on the detector. 

In order to separate the tunneled and reflected particles we have to take the limit t — > oo. Otherwise, the particles 
which tunnel after the time t would not contribute to the calculation. So we introduce operators 



F(x) = / dhD(x,ti). 
Jo 



N(x) 



dt\J (x, t\) 



(29a) 
(29b) 



From Eq. (g) it follows that the operator /t (oo, X) is equal to fp (X) + N (X). As long as the particle is initially 
before the barrier 



f T (X)p P (0) = p P (0)/ T (X) = 0. 



In the limit t — > oo we have 



Tun 



Tun 



Or) = 



1 



(x) = 



2{N(X) 
1 

2i(N(X) 



N (X) F (x) + F (x) N (X) ) , 
N(X),F(x) 



(30a) 
(30b) 



Let us define an "asymptotic time" as the integral of r(x,oo) over a wide region containing the barrier. Since the 
integral of tJ™ (x) is very small compared to that of r Tun (x) as we will see later, the asymptotic time is effectively 
the integral of T Tun (x) only. This allows us to identify r Tun (x) as "the density of the tunneling time" . 

In many cases for the simplification of mathematics it is common to write the integrals over time as the integrals 
from — oo to +oo. In our model we cannot without additional assumptions integrate in Eqs. ( p9| ) from — oo, because 
the negative time means the motion of the particle to the initial position. If some particle in the initial wave packet 
had negative momentum then in the limit t — > —oo it was behind the barrier and contributed to the tunneling time. 



VI. PROPERTIES OF THE TUNNELING TIME 



As it has been mentioned above, the question How much time does a tunneling particle spend under the barrier 
has no exact answer. We can determine only the time the tunneling particle spends in a large region, containing the 



barrier. In our model this time is expressed as an integral of quantity (30a) over the region. In order to determine 
the properties of this integral it is useful to determine properties of the integrand. 

To be able to expand the range of integration over time to — oo, it is necessary to have the initial wave packet far 
to the left from the points under the investigation and this wave packet must consist only of the waves moving in the 
positive direction. 

It is convenient to make calculations in the energy representation. Eigenfunctions of the Hamiltonian Hp are \E, a), 
where a = ±1. The sign '+' or '— ' corresponds to the positive or negative initial direction of the wave, respectively. 
Outside the barrier these eigenfunctions are: 

r | F , , v dfe ( exp (% PEX ) + r ^ exp (~& EX )) > x < °> (31a) 

' 1 JJ-t(E)exp(j-p E x),x>L, 
I^M—t (E) exp (-ip E x) , x < 0, 

<*!£,-> = < , ; " , „ ,n . (sib) 




(exp (-fax) - ^r* (E) exp (fax)) , x > L 
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where t (E) and r (E) are transmission and reflection amplitudes respectively, 



p E = V2ME, (32) 

the barrier is in the region between x = and x = L and M is the mass of the particle. These eigenfunctions are 
orthonormal, i.e., 

{E,a\E',a')=5 atCt ,5(E-E'). (33) 

The evolution operator is 

Up(t)=J2f o dE\E,a){E,a\exp^~Etj. 
The operator F (x) is given by the equation 

F(x) = J dhY, JJ dEdE'\E,a)(E,a\x)(x\E',a'}(E',a'\exp(^(E-E')ti\ 

where the integral over the time is 2irh5 (E — E') and, therefore, 

F (x) =2ttTiY^ J dE \E, a) (E, a\ x) (x \E, a') (E, a'\ . 

a,a' 

In analogous way 

N (x) = 2irhJ2 J dE \E, a) (E, a\ J (x) \E, a') (E, a'\ . 

a, a' 

We consider the initial wave packet consisting only of the waves moving in the positive direction. Then we have 

(a;)) = 2i:h J dE (\E, +) (E, +\ J (x) \E, +) (E, +\) , 
'P(x)N(X)) =4ir 2 h 2 Y, J dE (\E,+) (E,+\x)(x\E,a) (E,a\J (X)\E,+) (E,- 



From the condition X > L it follows 

'N 



(X))= JdE(\E,+)\t(E)\ 2 (E,+\). (34) 
For i<0wc obtain the following expressions for the quantities r Tun (x, t) and r^" (a;, t) 

(x, t) = / dE (\E, +> ±- \t (E)\ 2 (2 + r (E) exp (-2^) 



N(X)] 

evil ! ^ , 



+ r*{E)e^[2 % -p E x) ) (E,+\) , (35a) 



rS? (*>*) = i M v / ^ (l^ +> — I* (^)l 2 (E) cxp (-2^*) 
2 N(X)\J \ WE V V & / 



- r*(£) exp (2 -p B x ) )(E,+\). (35b) 



For x > L these expressions take the form 



* j dE {^ +> ^ |f (S)|2 ( 2 - ^ff exp { 2 ~n PEX 

^r(E)«p(-2i pj! x))^- 



T corr (""j 



2 (^7V (X) 



,/£' ( l-E, +) j- \t (E)\ 2 (jjj^r* (E) exp ^p E x 



(36a) 



(36b) 



t(E)^ E ^[- 2 n PEX )) {E '' 

We illustrate the obtained formulae for the ^-function barrier 

V (x) = Q5 (x) 

and for the rectangular barrier. The incident wave packet is Gaussian and it is localized far to the left from the 
barrier. 



2.0 



1.5 



Tun 
T 1.0 



0.5 



0.0' I 

-15 -10 -5 5 10 15 

X 

FIG. 2. "The asymptotic time density" for 5-function barrier with the parameter Q. = 2. The barrier is located at the point 
x = 0. Units are such that h = 1 and M — 1 and the average momentum of the Gaussian wave packet (p) = 1. In these units 
length and time are dimensionless. The width of the wave packet in the momentum space a = 0.001. 

In Fig. | and |, we see interferencedike oscillations near the barrier. Oscillations are not only in front of the barrier 
but also behind the barrier. When x is far from the barrier the "time density" tends to the value close to 1. This 
is in agreement with classical mechanics because in the chosen units the mean velocity of the particle is 1. In Fig 
|3|, another property of "tunneling time density" is seen: it is almost zero in the barrier region. This explains the 
Hartmann and Fletcher effect |2^,^] : for opaque barriers the effective tunneling velocity is very large. 

VII. THE REFLECTION TIME 

We can easily adapt our model for the reflection, too. For doing this, we should replace the tunneling-flag operator 
Jt by the reflection flag operator 

f R = 1 - It- (37) 
Replacing by fn in Eqs. ([50|) we obtain the equality 

(fa (t = co, X)) r Rcfl (x) = r Dw (x) - (/ T (t = oo, X)) r^ 1 " (x) . (38) 

We see that in our model the important condition 








FIG. 3. "The asymptotic time density" for rectangular barrier. The barrier is localized between the points x — and x = 5 
and the height of the barrier is Vo = 2. The used units and parameters of the initial wave packet are the same as in Fig. ^| 

r Dw _ Tr Tun + i?r Rofl (39) 

where T and R are transmission and reflection probabilities is satisfied automatically. 

If the wave packet consist only of the waves moving in positive direction, the density of dwell time is 

r Dw (x, t) = 2iTh J dE (\E, +) (E, +| x) (x \E, +) (E, +|) . (40) 

For i<0we have 

r Dw (x,t) =M J dE(\E,+)±- (l + \r(E)\ 2 +r(E)exp(~2^ PE x 

+ r*(E)exp(2^p E xj^(E,+\J (41) 

and for the reflection time we obtain the "time density" 

M 



T Hetl ^ 



i-(n(x)) 



I dE(\E,+)^-(2\r(E)f 



±(l + \r(E)\ 2 )r(E)exp(~2^ PE x\ + r* (E) exp I 2~p E r ) ) ,£.-|). (42) 



For x > L the density of the dwell time is 

r Dw (x,t)=M [ dE (\E, +) — \t (E)\ 2 (E 

J \ PE 

and the "density of the reflection time" may be expressed as 

W = f / dE (\E, + ) ± |t { E f { E) exp (2{p E x 

t* (E) 



t(E) 



(43) 



r (E) exp ( -2 l -p E x )) (E,+\) . (44) 



We illustrate the properties of the reflection time for the same barriers. The incident wave packet is Gaussian and 
it is localized far to the left from the barrier. 

In Fig. |and|, we also see the interference-like oscillations at both sides of the barrier. As far as for the rectangular 
barrier the "time density" is very small, the part behind the barrier is presented in Fig. pL Behind the barrier, the 
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FIG. 4. "Reflection time density" for the same conditions as in Fig. 




10 15 



FIG. 5. "Reflection time density" for the same conditions as in Fig. 



1.0x10 




-1.5x10 



FIG. 6. "Reflection time density" for rectangular barrier in the area behind the barrier. The parameters and the initial 
conditions are the same as in Fig. a 
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"time density" in certain places becomes negative. This is because the quantity r Refi (x) is not positive definite. 
Non-positivity is the direct consequence of non-commutativity of operators in Eqs. (|3^). There is nothing strange in 
negativity of r Rofl (x), because this quantity itself has no physical meaning. Only the integral over the large region 
has the meaning of time. When x is far to the left from the barrier the "time density" tends to the value close to 
2 and when x is far to the right from the barrier the "time density" tends to 0. This is in agreement with classical 
mechanics because in the chosen units the velocity of the particle is 1 and the reflected particle crosses the area before 
the barrier two times. 



VIII. THE ASYMPTOTIC TIME 



As it was mentioned above, we can determine only the time which the tunneling particle spends in a large r egion 
containing the barrier, i.e., the asymptotic time. In our model this time is expressed as an integral of quantity ( 30a ) 
over this region. We can do the integration explicitly. 

The continuity equation yields 

—D(x D ,t) + 7 —j(x D ,t) = 0. (45) 
at oxd 

The integration in Eq. ( [l9| ) can be performed by parts 

/■* - d r* 

J dtiD (x D ,ti) = tD (x D ,t) + — J txdtxJ (xn^tx) . 

If the density matrix pp (0) represents localized particle then lim^oo (^D (x, t) pp (0)^ = 0. The operator D (x, t) in 
all expressions under consideration is multiplied by pp (0). Therefore we can write an effective equality 







DC 



dtiD(x D ,h) = — J hdtxJ (x D ,ti) . (46) 
We introduce the operator 

I* oo 

f(x) = / hdhJix^h) . (47) 



We consider the asymptotic time, i.e., the time the particle spends between points x\ and xi when x\ — > — oo, 
x 2 — > +oo, 

t TUn (x 2 ,X 1 ) = / dxT TUn { X ). 

J X\ 

After the integration we have 

( X2 , X1 ) = t*™ (x 2 ) - t Tun (xx) (48) 

where 

(x) = -j± r (N (x) f(x)+f (x) N(x)). (49) 

2{n(x)) X 1 

If we assume that the initial wave packet is far to the left from the points under the investigation and consists only 
of the waves moving in the positive direction, then Eq. ( |48| ) may be simplified. 
In the energy representation 

f(x) = J t 1 dt 1 ^2jJdEdE'\E,a)(E,a\J(x)\E',a')(E',a'\exp(^(E-E')h 
Integral over time is equal to 2?7r?i 2 g§7<5 {E — E') and we obtain 
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f(x) = -ih'l-nh ^ J dE [ E, a) ^ • / J (r) E'.n') 



(E,a'\ 



E'=E 



+ (E,a\J(x)\E,a / }^(E,a'\^ , 
M^f^^iM^J^lE^^MJiXma) 



x (JL( E ,a\J(x)\E',+) 



E' = E 



+ {E,a\J(x)\E,+)—) (E,+\^). 



Substituting expressions for the matrix elements of the probability flux operator we obtain equation 



N(X)f(x) 



dE(*\E,+)t* (E)^t(E) (E,+\V) 

-Mx ( dE(V \E, +) — \t (E)\ 2 (E, +| *) 
J Pe 



+t nM- J dE(V \E, +) -^r* (E) t 2 (E) exp \2~Pex j ; E. -| * . 



When x — > +00, the last term vanishes and we have 



V(.V)T(,)) = / dE(*\E,+)t* {E)\^t (E)(E,+\V) 



^Mx j dE(^ \E, +) — It (E)\ 2 (E, + | x^+00. 

J PE 



(50) 



This expression is equal to < T (x) ) : 



N(X)T(x))^(T(x) / 
When the point with coordinate x is in front of the barrier, we obtain equality 



(51) 



XlX)T(.n) -ih j <IE(}V E. - )\l [ E)\ 2 ( J L ~ X 



M 



^. r (E)exp[~2p Ll ) + — ) (E.+\ m 



dE 



When \x\ is large the second term vanishes and we have 



N(X)T(x)) -^Mx { dE(y\E,+} — \t(E)\ 2 (E,+\y) 
J Pe 

+ J dE^\E,+)\t(E)\ 2 - i ^{E,+\^) 



(52) 



The imaginary part of expression (52) is not zero. This means, that for determination of the asymptotic time it is 
insufficient to integrate only in the region containing the barrier. For quasimonochromatic wave packets from Eqs. 
@, (H), ©, (|o|) and (|f) we obtain limits 



t 1 ^ n (x 2 ,x 1 ) ^t£ h + — M(x 2 -xi). 

PE 



(53a) 
(53b) 



where 



t^ = h—( & rgt(E)) 



(54) 
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X 

FIG. 7. The quantity r^"" (a;) for 5 function barrier with the parameters and initial conditions as in Fig. ^| The initial 
packet is shown with dashed line. 



is the phase time and 



t l ? = h^L{\n\t{E)\) (55) 



is the imaginary part of the complex time. 

In order to take the limit x — > — oo we have to perform more exact calculations. We cannot extend the range of 
the integration over the time to — oo, because this extension corresponds to the initial wave packet being infinitely far 
from the barrier. We can extend the range of the integration over the time to — oo only for calculation of N (X). For 
x < we obtain the following equality 

(jV(^)T(s)) = ^jf tdt(lUx,t)^h(x,t)-h(x,t)^IZ(x,t)) (56) 



where 

/■ (t i\ = I .//— . _ 

/Pe \n 



h{x,t) = I dE-^— \t (E)\ 2 exp (p E x - Et) I (E, + | (57) 



h(x,t) = J dE-^=(exp^p E x^j +r{E)exp(-^ PE x^cxp(^^E?j (£,+]*). (58) 

I\ (x, t) is equal to the wave function in the point x at the time moment t when the propagation is in the free space 
and the initial wave function in the energy representation is |t(-E)| (i?, +|\f r ). When t > and x — > — oo, then 
I\ (x, t) — > 0. That is why the initial wave packet contains only the waves moving in the positive direction. Therefore 

(^N (X)T (x)^ — > when x — > — oo. From this analysis it follows that the region in which the asymptotic time is 

determined has to contain not only the barrier but also the initial wave packet. 

In such a case from Eqs. (^8|) and (^) we obtain expression for the asymptotic time 

t Tun (x 2 ,x 1 ^-^)= f ^ f dE(*\E,+)t*(E)(—x 2 -ih-£-)t(E)(E,+\V). (59) 
N(X))J \Pe dEj 



From Eq. <pu) it follows that 



where T (X2) is defined as t he p robability flux integral (fl7|). Eqs. ([39j) and (|6^) give the same value for tunneling time 
as an approach in Refs. ]3l|,|^ 



{x%,Xi ► —00) = —j— r (T(x 2 )) (60) 

N(X)) 
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1.4 



1.2 - 




The integral of quantity t^" (x) over a large region is zero. We have seen that it is not enough to choose the 
region around the barrier — this region has to include also the initial wave packet location. We illustrate this fact by 
numerical calculations. 

The quantity t^™ (x) for (5-function barrier is represented in Fig. |7[ We see that r^™ (a;) is not equal to zero not 
only in the region around the barrier but also it is not zero in the location of the initial wave packet. For comparison, 
the quantity r Tun (x) for the same conditions is represented in Fig. |^. 

IX. CONCLUSION 

We have shown that it is impossible to determine the time a tunneling particle spends under the barrier, because 
the knowledge about the location of the particle is incompatible with the knowledge whether the particle will tunnel 
or not. This is because the corresponding operators, given by Eqs. (||) and ( |To| ) do not commute. However, it is 
possible to speak about the asymptotic time, i.e., the time the particle spends in a large region. 

In order to illustrate those facts, to obtain an expression of the asymptotic time and to investigate its behavior, 
we consider a procedure of time measurement, proposed by A. M. Steinberg |27[ ]. This procedure shows clearly the 
consequences of non-commutativity of the operators and the possibility of determination of the asymptotic time. 
Our model also reveals the Hartmann and Fletcher effect, i.e., for opaque barriers the effective velocity is very large, 
because the contribution of the barrier region to the time is almost zero. We cannot determine whether this velocity 
can be larger than c, because for this purpose one has to use a relativistic equation (e.g., the Dirac's equation). 

Due to non-commutativity of operators (Q) and ( |To| ) the outcome of measurements depends on particular detector 
even in an ideal case. This makes the measurement of the tunneling time difficult for opaque barriers, because the 
tunneling time is very short and the term depending on the detector increases linearly with the barrier width. This 
term vanishes when the time spent in a large region, including initial packet location, is measured. 
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